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Abstract 

According to a 2007 paper there was no evidence for a neutrino burst of two or more events in 
Super-Kamiokande (SK) during the entire period of data-taking from 1996 to 2005 from Andromeda 
or anywhere else. There is, however a scenario under which a detectable signal could have been 
missed given the search method employed by the analysis, and it would have been found using an 
alternate method. The alternate method depends on the hypothesis that two of the neutrino mass 
eigenstates have masses 4.0 eV and 21.4 eV which was inferred from an analysis of the SN 1987A 
data. Although one might argue that the hypothesis of such large neutrino masses is remote, there 
is a way they could be compatible with observed constraints on neutrino masses involving a third 
tachyonic {m? < 0) eigenstate, plus three sterile neutrinos. Given the importance of a positive 
supernova search result, and the ease of conducting it using existing SK data, there would seem to 
be little reason not to do it. 

PACS numbers: 13.15+g, 14.60Pq, 14.60St, 14.60Lni 
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I. INTRODUCTION 



To date the only supernova giving rise to an observed neutrino burst is SN 1987A for 
which 24 events were seen in 3 detectors in a burst lasting about 15 s llifl. Given the 
probable significant delay before future Mton neutrino detectors are completed, there is 
considerable interest in maximizing the efficiency for supernova detection using the best 
detectors currently available, especially the 50 kton Super-Kamiokande (SK). Two important 
ways of enhancing the sensitivity involve modifying SK in some way, or alternatively selecting 
an optimum energy threshold cut Efh to impose on the data that is based on the very different 
dependencies of the expected signal and background events 0, [sl- Of course, a supernova 
from a core collapse within our galaxy would be completely obvious in SK without needing 
to maximize its detection efficiency, and one well beyond Andromeda likely would not yield 
a single event in current detectors. This makes the possibility of seeing a supernova burst 
from Andromeda (expected signal 2 events in SK) one of great interest. Thus, given the 
spatial distribution of distances of potential sources, Andromeda plays a very special role in 
extragalactic supernova searches with SK. In fact, since the expected frequency of supernovae 
closer than a distance d does not increase significantly between the edge of our galaxy and 
Andromeda, and also between Andromeda (d = 778 kpc) and the next rich source at a 
distance of 2.7 Mpc, at which point the expected fiux would have dropped by a further 
factor of 10 by virtue of the inverse square law [g]. One may therefore conclude that if one 
should observe genuine signal of 2 - 3 events in SK the source is very likely to be Andromeda 
- a fact that could be checked based on directional information on those events. 

A 2007 paper analyzed data from SK for the period 1996 - 2005, and it reported no 
signal from Andromeda or anywhere else Q]. However, for the pre-2002 period SK had 
uncertainties in its energy calibration and it was not used in the 2007 analysis. While the 
pre-2002 SK could easily show a supernova in our galaxy, it was not sensitive enough to 
look for a supernova in Andromeda. Thus, the search for a signal from Andromeda, which 
required making an energy cut on the data in order to enhance the signal to background 
ratio, was based only on 3 years of data: 2005 - 2007. As of 2013 there now exists 7 more 
years of data to examine for SK, so it might be timely to redo the analysis to look for a 
low event rate signal from Andromeda. However, far more important than the additional 
7 years of data since 2005 is the possibility that the earlier analysis left open a way that a 
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detectable signal might have escaped notice. 



II. SEARCH METHOD USED IN 2007 ANALYSIS 

Given that the signal from SN 1987A lasted around 15 seconds, the search reported in 2007 
used a sliding 20 second time window, and it required the presence of at least N = 2 candidate 
events (a " cluster" ) in that window, nearly all of which would be expected to be background. 
In order to reduce background a threshold energy cut Et^ > 17MeV, was imposed which 
drastically reduced background to about one event per day, while simultaneously reducing 



any real supernova signal by about 33%, based on the expected supernova spectrum - 1 



lereby 



reducing the expected counts from a supernova in Andromeda from 2.0 to 1.4 js], [ij. The 
background was reduced further by considering the spatial positions of event vertices within 
the detector, since background events tend to cluster spatially, while real neutrino events 
would not. Therefore a cut imposed on the average spacing of the vertices for the clusters 
in any 20 second time window was found to reduce background almost to zero, while having 
a negligible effect on eliminating any real N event cluster In fact, the combined cuts 
on energy and vertex proximity left only 3 remaining candidate > 1 event clusters for 
the entire data-taking period. Unfortunately, each of those 3 clusters occurred within to 
6 seconds after blastings in the mine that housed SK, and they had to be excluded on that 
basis, thereby reducing the final number of candidates to zero. 

III. THE TWO MASS EIGENSTATE HYPOTHESIS 

Most analyses of SN 1987A assume that the 24 events can only set upper limits on the 



mass of the electron neutrino |8l-lll|. However, a recent analysis of SN 1987A neutrinos 



suggests that the (anti)electron neutrinos detected consist of two neutrino mass eigenstates 



having masses mi = 4.0 ± 0.5eV and m2 = 21.4 ± l.Oel^ 



explored much earlier 1^, 



3- 
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13| . using an idea that was 



Although the suggested masses seem preposterously large. 



the 2012-2013 papers suggest a way they could be consistent with much smaller empirical 
upper limits on the neutrino masses 1^, [l^. For example, the upper limit on the electron 
neutrino mass, and the sum of the flavor state masses from cosmology can all be addressed 
if the third mass eigenstate is tachyonic (m^ < 0) which would permit all three flavor states 
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to have masses very close to zero. It is also necessary to have 3 sterile neutrinos nearly 
degenerate with the three active ones to yield the well established Am^^^ and Am^j,; values. 

The claim of 4.0 eV and 21.4 eV mass eigenstates rests on the idea that most of the 24 
neutrinos observed were emitted during a time much less than the 15 second burst, possibly 
in as short a time as 1 second. Core collapse modellers differ on the time distribution in 
neutrino luminosity, and only some of them suggest a very large fraction in the first second. 
For example, while Totani et al. show 73% emitted in the first second and 82% in the first 
2 s (based on a numerical integration of their Fig. I),|l6| other sources suggest significantly 
less early luminosity [l3|. Thus, given the uncertainties, one cannot rule out the possibility 
that a large majority of the observed neutrinos from SN 1987A were emitted during an 
interval perhaps as short as 1 s. If one assumes near-simultaneous neutrino emissions then 
the neutrino arrival times at the detectors depend primarily on their mass m and energy E, 
based on relativistic kinematics according to: 

where t is the time a given neutrino is detected, to is the light travel time from the 
supernova, and T is some fixed time for all neutrinos associated with the supernova burst, 
i.e., the time t = T signifies when neutrinos of infinite energy E would have arrived. Note 
that equation 1 implies that on a plot of l/E"^ versus t all neutrinos having a fixed mass m 
lie on a straight line of slope M = 2/(m^to)- Thus, if neutrinos having two distinct masses 
were present among the 24 they would all lie on one of two straight lines having slopes 
Ml = 2/(rn'^tn) and M2 = 2/(m2to) that intersect at a point on the t axis. Remarkably, as 
Q, [3], this pattern is exactly what the data show (within the uncertainties in 
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shown in 

l/E"^) for every one of the 24 observed neutrinos. Finally, we note that the goodness of the 
two straight line fit (from which the claim was inferred) is not significantly impaired if there 
is a limited spread in emission times. In a Monte Carlo simulation in which 24 events are 
generated, assuming an emission time distribution and energy spectrum following Totanti 
et al.p^], and assuming the masses of individual neutrinos are either 4.0 eV or 21.4 eV, 
the agreement of the fit to two straight lines is usually just as good as for the actual data. 
Finally, it should be noted that the seemingly implausible claim of two such large mass 
eigenstates has recently received additional support from a sterile neutrino fit to the dark 
matter halos in the Milky Way and in galaxy clusters, using that same pair of masses jl^. 



IV. THE ALTERNATIVE SEARCH METHOD BASED ON 2 MASS EIGEN- 
STATES 



In the versus t plane the search window used in the 2007 search is a rectangle 

of height l/E"^ = 1/17'^MeV'^ and width 20 s which slides along the t axis. This search 
window is appropriate if we expect a neutrino burst from a source in Andromeda to have 
a similar width on arrival to that of SN 1987A, i.e., around 15 s. However, if the arrival 
times are primarily due to different travel times for neutrinos having different energies and 
masses, and if the masses are as large as 4.0 eV and 21.4 eV the arriving pulse will broaden 
the more remote the source, and individual neutrino arrival times will correlate with their 
energy according to Eq. (1). Since Andromeda is about 14 times as far as SN 1987A, a 15 
s burst from Andromeda would broaden to about 15 x 14 = 210s on reaching Earth. The 
alternate search method proposed here, therefore, replaces the fixed width rectangular search 
window by one consisting of two triangular regions - see Figure 1 - which bracket (within 
expected neutrino energy uncertainties) the two dashed lines having slopes corresponding to 
the masses 4.0 eV and 21.4 eV. Thus, the dashed line down the center of the larger triangle 
has a slope M = 2/{'m^to) where to = 2.54M?/ is the light travel time from Andromeda, 
and where m = 21AeV. Similarly, the dashed line bisecting the skinnier triangle has a slope 
given by the same equation but with m = A.OeV. The sides of the 2 triangles are based on a 
la ±16% uncertainty in E, and hence a ±32% uncertainty in 1/E'^ 20|. As the two-triangle 
search window slides smoothly along the t axis both triangular regions must of course move 
together to keep their bottom vertices in contact and on the t-axis. 

Note that unlike the fixed (20 s) width rectangular search window, the width of the two- 
triangle search window. At along the time axis is proportional to 1/E^f^, so that on a plot 
like Fig. 1, where the scale of the two axes is 1/E^^ and At, it is the rectangular search 
window that appears to change its width for two different values of Eth- For example, a choice 
Eth = 17MeV yields a width At = 63s whereas the choice Ef^ = lOMeV yields a width 
At = 182s for the two-triangle window. It is clear that the proposed two-triangle search 
window might indeed produce candidate events that would be missed by the standard search 
window - especially when a threshold Eth = lOMeV were used instead of Efh = l7MeV and 
the actual area of the search window becomes nine times larger. As an example, consider the 
3 simulated events shown in Figure 1 as black circles and assume that their locations in time 
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are appropriate to the choice At = 182s. They clearly are all inside the two-triangle search 
window with Eth = lOMeV, but only one of the 3 would be picked up as a candidate either 
for a rectangular search window or for a two-triangle search window with Eth = 17MeV. 
Hence no > 2 event cluster would register in this case using either of these latter two 
search windows - the standard rectangular one or the two-triangle one with too high an 
energy threshold. 

One might argue justifiably that using an energy threshold as low as Eth = lOMeV could 
generate a number of false signals. However, while Eth = 17MeV may be the optimum 
choice in terms of enhancing the signal to background rates, it may not be an optimum 
choice if very few real events are expected as in the case of a signal from Andromeda. If we 
expect only two events in SK from a supernova in Andromeda, then according to Poisson 
statistics there is a 27% chance of getting 2 events and a 33% chance of getting 3 or more 
events. In contrast, if a cut Eth = l7MeV is applied the expected number of events would 
drop to 1.4 instead of 2.0, and the previous numbers become: a 24% chance of getting 
2 events and a 17% chance of getting 3 or more events. This means that the chances of 
finding a very persuasive signal (A^ > 2 events) essentially is reduced in half by virtue of the 
Eth = 17MeV cut. Furthermore, given that zero events were seen in the original search, the 
possibility of a few background events masquerading as a signal is less worrisome when one 
might easily eliminate them if they do not point back to Andromeda. The authors of the 
original search did in fact try lowering Eth after finding zero events with the original choice, 
and they still found no signal. The main argument being advanced in the present paper, 
however, is not the that a signal could have been missed due to too low an Eth = 17MeV 
cut, but rather due to the rectangular shape of the search window. 



V. WHAT IF A SIGNAL IS SEEN USING THE ALTERNATIVE METHOD? 

The frequency of supernovae in a galaxy is very uncertain, since perhaps as many as 
80 — 90% are hidden by dust [3l • Andromeda is a galaxy of comparable size to the Milky 
Way, and the expected frequency of supernovae is also comparable - perhaps 5 per century, 
or one every 20 years - which suggests there could be a 50% chance that a supernova 
occurred during observing period without any visible indication of it. If no signal is found 
by examining the 2002-2012 record for SK using the alternate search method, it could 
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simply mean that there were no supernovae in Andromeda during the 10 years since 2002. 
If 2 to 3 events are observed in SK, one would need to see how well they point back to 
Andromeda. At an energy of lOMeV the directional uncertainty of arriving neutrinos is 



about 12° 



20|. A directional uncertainty in the source of about 12° amounts to about 



1% of all possible solid angles, so the chances of finding a cluster of N events all pointing 
back to Andromeda within 12° would be p = 0.01^ which would be highly convincing for 
> 1. Were such a result found and the events not picked up within the earlier 20 s wide 
rectangular search window, it would constitute evidence for a pulse broadening with distance, 
in conformity with the existence of 4.0 eV and 21.4 eV neutrinos. Hence, the significance 
of such a positive result would be far greater than being the first extragalactic observation 
of a supernova burst of neutrinos, since it would also give support to the hypothesis of two 
identified mass eigenstates, and all that follows from that hypothesis, including the existence 
of tachyonic neutrinos. Such a hypothetical positive search result would constitute a third 
piece of evidence in support of the two claimed mass eigenstates in addition to the SN 1987A 



analysis 12|, ll3| and the sterile neutrino dark matter fits. 19|. 
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FIG. 1: On a plot of versus t the figure sliows tlie searcli windows used in the original 

2007 analysis (solid and dashed rectangles), and the "two-triangle search window" that is based 
on the assumption of 4.0 eV and 21.4 eV mass eigenstates in SN 1987A. Note that as the 
two-triangle search window slides along the t axis, it would contain the 3 simulated events 
for the position shown. The value of the total horizontal width along the time axis At of the 
two-triangle search window is proportional to ^/E^/^ and it is 182 s for Eth = WMeV, and 62s 
for Efh = 17MeV. When Ef^ = 10(17)Mey, the 20 s wide rectangular search window would 
be the skinny (fat) rectangle. Hence, we see that both the fat and the skinny rectangles would 
contain only one event at a time (as it slides along t), and it would miss the hypothetical 
signal. Moreover, the two-triangle search window with the choice Ef^ = 17MeV would also miss 
the signal because in this case the 3 events would not be as depicted in the figure, but would 
be spread out in time three-fold from their depiction there, when the time axis is in units of seconds. 
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